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Abstract 

The scattering properties of the non-hnear 0(3) model in (2+l)-D, 
modified by the addition of both a potential-hke term and a Skyrme- 
hke term, are considered. Most of the work is numerical. The skyrmion- 
scattering is found to be quasi-elastic, the skyrmions' energy density 
profiles remaining unscathed after collisions. In low-energy processes 
the skyrmions exhibit back-scattering, while at larger energies they 
scatter at right angles. These results confirm those obtained in pre- 
vious investigations, in which a similar problem was studied for a 
different choice of the potential-like term. 



1 Introduction 

In the past few years, a-models in low dimensions have become an increas- 
ingly important area of research, often arising as approximate models in the 
contexts of both particle and solid state physics. They have been used in the 
construction of high- Tc superconductivity and the quantum Hall effect; in 
two Euclidean dimensions, they appear to be the low-dimensional analogues 
of four- dimensional Yang-Mills theories. Moreover, they are examples of the 
harmonic maps studied by differential geometers and, as such, are interesting 
in themselves. But only very special a-models in (2+l)-D are integrable [Q], 
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and the physically relevant Lorentz-invariant models are not amongst them; 
in these cases, recourse to numerical evolution must be made. 

The simplest Lorentz-invariant model in (2-|-l)-D is the 0(3) model, 
which involves three real scalar fields, (f){x^)={(f)a{x^) , a=l,2,3}, with the 
constraint that (f) lies on the unit sphere 32^^ : 

0.0=1. (1) 

Subject to this constraint, the Lagrangian density and the corresponding 
equations of motion are 

= \{d,cf>)W^), (2) 

d^ct> = [-{dtd^f + {d,(t>f + {dy<t>f\<t> + dlct> + (3) 

Note that we are concerned with the model in (2-|-l)-D: x'* = {x^,x^,x^) — 
{t, X, y), with the speed of light set equal to unity. 

An alternative and convenient formulation of the model is in terms of one 
independent complex field, W, related to the fields 4>a via 

^ - (4) 

In this formulation, the Lagrangian density and the corresponding equations 
of motion read (asterisk denotes complex conjugation) 



and 



The problem is completely specified by giving the boundary conditions. As 
usual we take 

lim 4>{r, 9, t) = Mt), (7) 

where (f>o{t) is independent of the polar angle 6. In (2+l)-D this condition 
ensures a finite potential energy, whereas in two Euclidean dimensions, i.e., 
when (f) is independent of time, it leads to the finiteness of the action, which 
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is precisely the requirement for quantization in terms of path integrals. As 
shown by several people 0, any rational function W{z) or W{z*), where 
z = X + iy, is a static solution of Eq. (^). These are the instantons of the 
model, and can be regarded as static solitons of the same model in (2+l)-D. 
The simplest one-soliton solution, W = Xz {X is a free parameter determining 
the size of the soliton) has been numerically studied by Leese et.al When 
viewed as an evolving structure in (2+l)-D, the soliton has been found to 
be unstable. Any small perturbation, either explicit or introduced by the 
discretization procedure, changes its size. This instability is associated with 
the conformal invariance of the 0(3) Lagrangian in two dimensions. 

The 0(3) solitons, however, can be stabilized through a judicious intro- 
duction of a scale into the model, thereby breaking its conformal invariance. 
This has been done by Leese et.al. p and, using a more general field, by 
ourselves 0. 

In the present article we present the results of soliton-like scattering in 
(2+l)-D obtained by applying the methods of [^] to a more general skyrmion 
field. In the next section we present our skyrme model and give a brief 
account of our previous paper where the same model was considered for 
the case of zero-speed systems. After explaining the numerical procedure in 
section 3, we pass on to discuss two-skyrmion scattering in section 4. The 
closing section contains our conclusions. 



2 Skyrme model in (2+1) dimensions 

Using the VT-formulation, our Skyrme model is defined by the Lorentz- 
invariant Lagragian density 

C = 

no XdtW*dyW - dtWdyW*f + {dtW*d,W - dtWd,W*f 

" ^^'^ {i + WW 

{d^W*dyW -d^WdyW*f 

(i + |iyp)4 ^ 
\w - x\^ 

- 4^2 , „,.,L , (8) 



where is given by Eq. (^) and ^1,^2 are real parameters with dimensions 
of length squared and inverse length squared, respectively; they introduce a 
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scale into the model, which is no longer conformal invariant. If the size of 
the solitons is appropriately chosen, it is energetically unfavourable for the 
solitons to change it. The ^i-term is the (2+l)-D analogue of the Skyrme 
term, whereas the ^2-term is a potential-like one. Unlike the former, the 
latter term is highly nonunique 0. 

The field equation corresponding to the above Lagrangian can be cast 
into the form 



"yy ^ 



■[2WlWtW, + 2W*WtWy - 2W*WMy 



(^X _)_ ||y |2^2 L-' ' ti: " i-" ^- ' • ty'" y ^""'xy' 

+ w:m^ - w^) + w;y{wl - w^) - w;,{wl + 

+ W,,{\Wt\^ - \Wy\^) + Wyy{\Wt\^ " \W ^\^) + W u{\W + \W y\^) 

+ w,y{w:Wy + w;w^) - Wtxiw^w, + w:wt) - Wty{w;Wy + w;wt) 

2W 

+ YV^2^^^tWy - w;w,f + {wiwt - wiw.f - {w;Wy - w;w,y 



+ n ...m . (9) 



16^21 


\W - 


- A| 


2 


(1 + 


\W\ 


2^3 



where the notation Wx = dxW, Wxx = d^W, etc., has been used. It is 
straightforward to check that 

W = (10) 

z — b 

is a static solution of Eq. (y), provided the following relation holds: 



A^^. (11) 

a — b 

This is the familiar one-instanton of the non-linear 0(3) model, but now 
A, which characterizes its size, is no longer a free parameter: It is fixed by 
Eq. (|TTD . The instanton, with its size thus fixed, is usually referred to as 
a 'skyrmion'. One can readily check that the maximum of its total energy 
density is given by 

(\\\2 I 1 \2 

Ema. = e{l + e,e), e = 8||^ (12) 

|A(a-6)|2 
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For the parameter values given below in section 3, Eq. (|1^) yields Emax = 
129.3, the 'canonical size'. The numerically-obtained Emax is very near this 
value 

(See Figure 1). The distance from Emax to the centre of the lattice is given 
by 

(a|A|2 + 6)/(|A| + l). 

To study scattering processes we take the field 

W = A^^. (13) 

which describes two instantons. However, as there are interaction forces 
between the instantons, Eq. ( [TB| ) is not a static solution of the equations of 
motion. 

In our previous paper where we studied only the case of skyrmions 
started off at rest, we numerically evolved Eqs. (|10D and (|13]) and found two 
basic results: 

• The field (|T0[) remains almost perfectly static, its total energy density 
being practically unaltered as time elapses (see Figure 1); 



• The field (|T3[) shows two skyrmions shaking off some kinetic energy, 
thus adjusting themselves to their canonical sizes. This is in accordance 
to expectation, as this field is not a solution of the equations of motion. 
Then the skyrmions slowly move away from each other, unveiling the 
presence of a repulsive force between them. 

These observations confirm the results obtained in 0, where the skyrmion 
of the model was just \z. This configuration posseses a total energy density 
whose maximum is positioned at the centre of the grid i^Zmax = 0), the ana- 
logue of 
Eq. ([1|) being Emax = 8X^(1 + S^aA^), A = {/ 261/62. With regards to 



the scattering, reference considered fields of the form W = X{z ~ a){z + b) 
which, as their cousins (|13D, resemble two interacting instantons in an ap- 
proximate manner. 
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3 Numerical procedure 



Our simulations employed the fourth-order Runge-Kutta method, and ap- 
proximated the spatial derivatives by finite differences; the Laplacian was 
evaluated using the nine-point formula. Worthy of note is the fact that the 
finite-difference expressions for the derivatives of the fields \z and \{z — 
a){z + 6), used in 0, are exact. This is no longer true for our more gen- 
eral model which, in this sense, is more 'perturbed'. Fortunately, this factor 
turned out to have no significant effect on the qualitative properties of the 
discretized version of our model. 

All computations were performed on the workstations at Durham, on a 
fixed 201x201 lattice with spatial and time steps 5a;=5?/=0.02 and 5t=0.005. 
Every few iterations we rescaled the fields (p (f)/ \/ (p.cf) . We also included 
along the boundary a narrow strip to absorb the various radiation waves, thus 
reducing their effects on the skyrmions via the reflections from the boundary. 
As time elapses, this absorption of radiation manifests itself through a small 
decrease of the total energy, which gradually stabilizes as the radiation waves 
are gradually absorbed. 

We choose the parameters to have the values: ^1=0.015006250, ^2=0.1250, 
a = c = 0.75 and b = d = 0.05 which, according to Eq. (PT]), set A=l. 



4 Skyrmion scattering 

Our simulations look at the initial conflguration 

.-0.75^ + 0.75 

z - 0.05 z + 0.05' ^ ' 

and evolve it for different initial velocities. First, let us consider head-on 
collisions. There is always an initial burst of radiation as the skyrmions 
regain their canonical size. At small speeds, the skyrmions approach each 
other, but the repulsive force between them results in their motion being 
reversed. In Figure 2 we present some pictures of the total energy density 
for skyrmions with initial speed equal to 0.2; the corresponding contour plots 
are shown in Figure 3. 

A qualitatively similar behaviour is observed for speeds up to approxi- 
mately 0.3, for which the skyrmions acquire enough kinetic energy to over- 
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come their mutual repulsion; during their collision they form a rather com- 
plicated state and then re-emerge at 90° to the original direction of motion. 
The emerging skyrmions are initially shrinking but, after they have travelled 
some distance, they expand once more. This final state is achieved after some 
oscillations of the energy density. Figure 4 and Figure 5 show several total 
energy density pictures and contour plots for this 90°-scattering, whereas 
Figure 6 exhibits how the amplitude of the total energy density for the above 
cases varies in time. 

In the case of non-zero impact parameter the results are very much as 
expected. With an impact parameter small enough to prevent the skyrmions 
from getting too quickly to the boundary, they scatter almost backwards 
along their original trajectories or at 90°, depending on the initial speed. 
Also, the larger the impact parameter, the smaller the scattering angle. A 
typical case is shown in Figure 7. 

On comparing these results with those obtained for the configuration 
W = X{z — a){z + b) we find there is no qualitative difference, and hence 
the results obtained in are borne out by the general case studied in the 
present work. 



7 



Emax=128.5 





Emax=1 14.4 




Figure 2: Total energy density pictures corresponding to the initial velocity 
V = (0.2,0.0). 
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Emax=134.9 




Figure 2: Continued. 
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Figure 4: Scattering at 90 degrees for the case v — (0.3, 0.0). 
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Emax=128.2 
t=2.5 




13 







2 







2 



t=2.5 





Figure 5: Contour plots for the scattering of Figure 4. 
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Figure 6: Maximum of total energy density vs. time. 
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t=2.5 



t=5 













Figure 7: Non-zero impact parameter collision. The initial velocity is ^; = 
(0.3,0.01). 
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5 Conclusions 



We have performed a numerical study of the scattering properties of the 
general two-instanton (-skyrmion) field configuration of our version of the 
Skyrme model in (2+l)-D, confirming the results previously found in other 
version of the model where, unlike the case studied in this paper, the finite- 
difference expressions for the derivatives of the fields enjoyed the numerically 
convenient feature of being exact. This factor, however, did not affect the 
qualitative scattering properties of the model under study. 

Although our field configurations resemble two instantons in an approx- 
imate manner -the model is not integrable-, they exhibit a clear soliton-like 
behaviour. All radiation effects are small and the skyrmions' total energy 
density profiles are preserved during the scattering process. They may get 
distorted, but always recover when the distance between the skyrmions be- 
comes large enough. Their interaction is of a repulsive nature, at least at large 
distances, the interaction being more difficult to asses when the skyrmions 
are close together. For head-on coUisions there is a velocity below which 
the skyrmions bounce back, and above which they scatter at ninety degrees. 
There is a resemblance both with the properties of kinks in the 0^ model, 
which has a critical velocity, and with monopole scattering at 90°. 
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